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^ ^ A good generalization of the Euclidean dimension to disordered systems and non crystalline 

, structures is commonly required to be related to large scale geometry and it is expected to be 

' independent of local geometrical modifications. The spectral dimension, defined according to the 

low frequency density of vibrational states, appears to be the best candidate as far as dynamical 
and thermodynamical properties are concerned. In this letter we give the rigorous analytical proof 
of its independence of finite scale geometry. We show that the spectral dimension is invariant 
under local rescaling of couplings and under addition of finite range couplings, or infinite range 
• couplings decaying faster then a characteristic power law. We also prove that it is left unchanged 

by coarse graining transformations, which are the generalization to graphs and networks of the 
usual decimation on regular structures. A quite important consequence of all these properties is 
the possibility of dealing with simplified geometrical models with nearest-neighbors interactions to 
study the critical behavior of systems with geometrical disorder. 
. On translationally invariant structures such as crystals, many physical phenomena fundamentally depend on the 
c/2 Euclidean dimension d, regardless of the particular form of elementary cells. This parameter appears to resume all 
^ ,. relevant informations about large scale geometry and therefore it affects the dynamics at low frequencies as well as 
■ the thermodynamics at low temperatures and near critical points. In fact this relation is one of the deepest aspects 
^ , of the physics of systems with many degrees of freedom. Even if it has not been proved yet in general, it is commonly 
' , ■ accepted and it is experimentally found that this dimensional dependence is not modified when translational invariance 
^ is lost as a consequence of inhomogeneities or disorder at finite scales. This invariance under couplings disordering and 
Q local changes in geometry, which we call geometrical universality, is perhaps the most important feature of dimension 
^ in statistical physics. Therefore it is quite reasonable to require that it is preserved in any meaningful generalization 
, of the Euclidean dimension. 
I ' The spectral dimension d, firstly introduced by Alexander and Orbach [Q to characterize the low frequency vibra- 
^ J tional spectrum of fractals, is considered for many reasons the right generalization of d to non-crystalline structures 
QQ ' such as polymers, glasses, fractals, amorphous materials and so on. In fact d exactly replaces d in most laws where di- 
' mensional dependence explicitly appears: the spectrum of harmonic oscillations, the average autocorrelation function 
I of random walks, the critical exponents of the spherical model |^], the low temperature specific heat, the generalized 
Mermin- Wagner theorem the infrared singularities of the Gaussian model and many others. In addition d can be 
easily measured by well known experimental techniques such as e.g. neutron scattering. These remarkable properties 
suggest extending to non-crystalline structures the same picture holding for crystals. In other words, we expect that, 
even on non translationally invariant networks, the fundamental physical properties depend only on large scale geom- 
' etry and that all this dependence is encoded in d. This would mean that all we need to deal with geometrical disorder 
' is a non-integer parameter generalizing the usual dimension and it also would mean that geometrical disorder does 
I not imply complex behavior. Clearly, such an appealing picture strongly depends on universal properties of d which 

are not obvious neither proved yet. The rigorous proof of these properties is the subject of this letter. 
^ , After introducing the basic concept and notations concerning J, we will prove its geometrical universality following 
three major steps, modifying the networks by acting on the strength of links, on the number of links and on the number 
of points. First, d is proved to be invariant under a bounded local rescaling of couplings, which therefore does not 
alter the underlying geometry. Then we will prove its invariance under the addition of finite range and infinite range 
couplings up to a characteristic power law decay, changing in this way the geometrical structure under consideration 
by modifying an infinite number of links. Finally we will consider a general coarse graining procedure, which changes 
[ the number of degrees of freedom of the network while leaving d invariant. These general transformations allow to 
group very different networks under the same geometrical universality class, characterized by a given value of d. Their 
significance and applications will be discussed in details in the conclusions. 
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The spectral dimension of an infinite connected graph G can be defined by studying the low frequency behavior 
of its vibrational spectrum. A graph is a set of points with a topology described by its adjacency matrix , whose 
elements are equal to 1 if the points i and j are nearest neighbors and equal to in all other cases. 

This mathematical structure can be used as a geometrical model for an oscillating network with masses m at 
each point, connected along the links by springs with elastic constant k. Introducing the positive definite Laplacian 
operator Lij = Zi5ij — Aij, Zi — 'Y^- Aij being the coordination number of site z, the eigenvalue equations for the 
vibrational normal modes read 



kY^L 



ijjLj = uP'm Xi (1) 



where u) is the frequency and Xi the displacement from equilibrium position at site i. The spectral dimension d is then 
defined by the asymptotic behavior of the density of modes with frequency lo, Pui^) ~ w''"^ for w — > 0. By (^), the 
spectral dimension is related to the density pi{l) of eigenstates of the Laplacian with eigenvalue I, whose behavior is 
given by pi{l) - l"^/^-^ for / -> 0. This relation gives the basis to describe large scale properties of statistical models 
in terms of d. In particular, we shall consider the Gaussian model, which provides the most effective tools to deal 
with universal properties of d. 

The Gaussian model on G is defined by the Hamiltonian: 



and its specific free energy / is given by 



/ = - lim ^ F^ - lim ^ logZ (3) 

where Z is the partition function calculated according to the Boltzmann weight exp{—H). Its correlation functions 
are given by 

<<P.cl,,>^{L + m%' (4) 

and from the matrix <I>y =< (j>i(l)j > is positive definite. The spectral dimension is related to the singular part of 
/ for ^ by: 

Stng if) ^ m^. (5) 

Using this definition of d and the above mentioned properties, a "weak" universality of d has been proved, namely its 
independence of a local bounded rescaling of masses [Q, showing that is left invariant replacing by mf = aim^ 
in with K^^ < < K for some K > 1. 

Now our first step is the extension of such invariance to a local bounded rescaling of couplings. Such "stronger" 
universality will provide the mathematical basis for all further points. First of all, we notice that a global rescaling of 
all couplings by a given constant c does not affect d, since it is equivalent to a global mass and correlation functions 
rescaling. Then we prove the invariance of (|^) under the rescaling Aij —> Jij = CijAij, with < kij < K, 

corresponding to the generalized Laplacian L'^^ = z[5ij — Jij, z'i — Jy. The strategy consists in taking successive 
derivatives of (||) with respect to rr? up to a divergence for ^ 0, in proving the monotonicity of such derivatives 
with respect to a generic local increasing of the couplings Jy and then exploiting their boundedness between the two 
limit cases of global rescaling = and Jy = K . The key point is the monotonicity relation 

holding for a^i > and every n>\. This can be proved by noticing that 

d ( d 



dJki V dm^ 



F = (<i>")fefe + - 2($")fe, (7) 



and that the right hand side of is the value of a positive definite quadratic form computed on two coinciding 
vectors v, namely Y^^. Vi{^'^)ijVj, with Vi = Sik ~ Su. 
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This first theorem (which we will call strong universality theorem or SUT) proves that d depends only on topology 
and is not affected by possible disorder in interactions strengths and, together with the previous result [Q concerning 
rescaling of masses, shows that the low frequency vibrational spectrum is only related to geometry. In the following, we 
will see that indeed only large scale geometry comes into play. The second and third steps deal with transformations 
affecting the topology itself, implying a change in the number of links and points. It is clear that the changes we are 
considering must involve sets of links or points having non vanishing density in the whole graph in the thermodynamic 
limit. Indeed, changing a finite number of entry in the Laplacian, or an infinite one with zero density, corresponds to 
a compact perturbation and therefore does not modify its continuous spectrum |5| . 

The simplest transformation on links is the addiction of bounded next to nearest neighbor interactions (NNI). Their 
exact values are clearly irrelevant due to SUT. A possible system G' with such additional interaction is described by 
the generalized Laplacian L'^^ — cLij — {L'^)ij, with c > 2Zjnaxj which is second degree polynomial in lacking of the 
constant term. Its eigenvalues I' are related to the eigenvalues I oi L hy I' — kl ~ P and they become asymptotically 
proportional for / 0. This implies the asymptotic coincidence of the respective spectral density and, therefore, the 
invariance of d. Now, due to SUT, a generic graph obtained from G by adding some (not necessarily all!) NNI has 
the same spectral dimension too. This proof is easily extended to the addition of interactions up to an arbitrary finite 
distance q. Indeed, iterating n times the same NNI transformation on G', one obtains a system with interactions 
to all distances up to 2"+^. Choosing a suitable n, the invariance of d follows from SUT. Now, consider the inverse 
transformation, consisting in acting on G cutting bonds which are limited range couplings with respect to the topology 
of the resulting graph G". We shall call it bond cutting transformation (BCT). From its definition, it follows that BCT 
preserve d. This property will be very useful in the following. 

An analogous strategy can be used to show the invariance of d under the introduction of a long range interaction, 
with power law or exponential decreasing in the graph chemical distance r. The fundamental observation is that, 
when the asymptotic behavior of the eigenvalues V of the generalized long range Laplacian L' is given by /' ^ cl + 1", 
with a > 1, the new interaction does not modify d, since the above exploited proportionality still holds. 

Now, if the long range interaction decreases as r~'^ , where ^ > dc {d^ being connectivity dimension of G) as required 
by local energy finiteness, it can be shown, expanding L' on the basis of the eigenvectors of L, that a = d{'y — dc)/2dc. 
Therefore J is left unchanged by power law interactions with 7 > dc + 2^ and, from SUT, by any exponential decaying 
interaction. Notice that on Euclidean lattices, where dc — d ~ d, this reproduces the well known result 7 > d + 2. 
The addiction of a power law interaction with a slower decay always alter d. From the previous considerations, d 
increases according to d' = 2dc/ (7 — dc). 

The class of transformations considered up to now does not modify the number of points of G. 

The last kind of graph modifications we will analyze involves a non trivial change in the number of degrees of 
freedom. Actually they include one of the most important class of transformations used in statistical mechanics, 
namely decimations, which are the basic tool of renormalization group on lattices and fractals. 

We call these transformations topological rescalings, since they generate local changes in the number of points. The 
most general topological rescaling can be realized through two independent steps. The first one is the partition and 
consists in dividing the graph G in an infinite family of connected subgraphs Ga, with uniformly bounded number 
of points. The second one is the substitution and consists in generating a new graph G' by replacing some or all Ga 
by a different (connected) graph Sa, whose number of points ranges from 1 to a fixed Nmax, and by adding links 
connecting different Sa in such a way that two generic Sa and Sfs are connected by some links if and only if Ga 
and G/3 were. The simplest topological rescaling occurs when every Sa is composed by just one point. In this case 
the resulting graph Gm is called the minimal structure of the partition {Ga}- The invariance of d under topological 
rescaling can be reduced to proving that G and the minimal structures of its partitions have the same d. Indeed, from 
this property it easily follows that two graphs having the same minimal structure also have the same d and it can be 
proven that G and G' have the same minimal structure. We will now show that d is left unchanged reducing G to its 
minimal structure. 

For every subgraph Gq let us define the boundary of Ga , dGa , as the set of points of Ga which are connected to at 
least one point oi Gf3, (3 ^ a. We will call Ma the number of points of dGa- Let us consider a new graph G' obtained 
by replacing each subgraph Gq with the complete graph of its boundary points, Km^ (the complete graph Kn is the 
n-points graph with each point connected to all the other ones). It can be shown B that, considering the Gaussian 
model (|2|) on G and the Gaussian model on G' with suitable bounded masses m'^, defined on the points of the dGa 's 
and bounded couplings J^ij, between the points of Km^, i one has for the free energies defined on the two graphs G 
and G': 

/g = /g' (8) 

The geometry of the new graph G' can be further simplified applying the previous transformations. First, since we 
are interested in the singular part of /, the exact values of the bounded masses and couplings between the points of 
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Km^ are irrelevant due to SUT. 

Now consider a single Km^, together with the Km^i which are connected to it. Then for every Km^ select one point 
iaifi) among the Ma{(3) points connected to and cut the links between all the other Ma{(3) — 1 points and KMfj 
using BCT. 

Applying again the bonds cutting transformation, one can reduce the ia{P) points to a linear chain, observing that 
all the other links can be considered as couplings up to k-th neighbors with respect to nearest-neighbors interactions 
between points of the linear chain. 

The linear chains can now be shrinked to one point pa by first applying the previously mentioned Gaussian model 
reduction [^, leading to a new relation for the free energies as in (^), and then cutting the spurious links between 
external points added with this operation using BCT. 

Therefore, one is left with a single point Pa for each Gq, connected to the pfj corresponding to those subgraphs 
G/3 whose points where connected to Ga in G. This graph Gm precisely corresponds to the minimal structure of the 
initial partition of G and it has the same d as G. 

The three very general classes of geometrical transformations we have been considering can be applied in all possible 
sequences to a graph, leading to an overall transformation on coupling strength, number of links and degrees of freedom 
which does not change its spectral dimension d. We will call such a transformation an isospectrality. Interestingly 
enough, isospectralities include a well known class of transformations on graphs studied in graph theory, namely the 
rough- isometrics [Q. 

Indeed, isospectralities include most part of currently used transformations. As an example, the usual decimation 
procedure on fractals is a topological rescaling. In particular, for all exactly decimable fractals (such as e.g. Sierpinski 
gaskets and T-fractals the minimal structure of the graph coincides with the graph itself. Again, an isospectrality 
relates the usual two dimensional square lattice, the hexagonal lattice and the triangular lattice, which therefore all 
have dimension 2. In other words, isospectralities are the theoretical formalization of the intuitive idea of invariance 
with respect to bounded scale perturbations and disorder and the isospectrality classes, defined as the classes of 
graphs related by such transformations, are the practical realization of the apparently abstract concept of non integer 
dimension. Now, since most dynamical and thermodynamical properties of generic discrete structures depend only on 
d, isospectralities provide a very powerful tool to reduce a very complicated geometrical structure to the simplest one 
having the same d. The latter turns out to be much simpler to study and still presents the same universal properties. 
Moreover, not only an isospectrality can be used to reduce and simplify structures and problems. It can also be 
applied, with the opposite aim, to build complicated structures with controlled dynamical and thermodynamical 
properties, starting from simple deterministic geometrical models. This is the point of view of spectral dimension 
engineering, providing a very interesting field of possibilities to polymer physicists and material scientists dealing 
with non-crystalline materials. 
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Figure caption 

Fig.l: Example of isospectral structures obtained applying isospectral transformations (without long range cou- 
plings) to the T- fractal (a) and to the square lattice (b). 
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